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Laboratory tests conducted to determine the effects of wave groups on rubble 
mound breakwater stability are analyzed. A literature review outlining the controversy 
regarding whether or not wave groups effect breakwater stability is given. Resolution of 
this controversy is proposed using two independent parameters that characterize wave 
groups. The two parameters are an envelope exceedance coefficient a and spectral 
shape y. The incident time series must be known to compute a, and an algorithm is 
given that will resolve incident and reflected wave time series from closely spaced wave 
gauges. The digital to analog simulation of laboratory waves with specified wave group 
characteristics is reviewed. It is recommended that further studies be conducted to 
parameterize the envelope exceedance coefficient a with a non-biased groupiness factor. 
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INTRODUCTION 



The design of rubble mound breakwaters is one of the oldest and most common 
applications of coastal engineering methodologies. Design of rubble mound structures is 
most commonly based on the empirical methods presented in the Shore Protection 
Manual (SPM 1984). Due to some failures of structures designed with the SPM 
methodology, an effort to improve the SPM design methodology is undertaken. The 
design methodology presented in the SPM is based on empirical formulas developed 
from experiments with monochromatic waves. Most of the more recently proposed 
design criteria are based on laboratory experiments using irregular waves. The effects of 
spectral shape and wave grouping will be targeted. As a proposed improvement, factors 
such as permeability, wave period, and storm duration may have secondary effects on 
rubble mound breakwater stability. 

It has been well-documented that waves in random seas tend to form groups not 
accounted for by the purely Gaussian random wave model. Many of the coastal 
engineering problems associated with wave groups have been documented by Medina 
and Hudspeth (1990). In the last two decades, controversy has arisen as to whether 
spectra composed of grouped waves cause different levels of damage to rubble mound 
structures than spectra of ungrouped waves. Medina et al. (1990) have found that most 
of the controversy may be resolved if two independent spectral related parameters are 
chosen to characterize the wave groups (Mase and Iwagaki 1986). One parameter 
characterizes the length of wave groups and the other characterizes the magnitude of the 
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energy flux in a group of waves. 

Intermediate scale laboratory studies to determine if wave groups affect rubble 
mound breakwater stability were performed in the 2-dimensional wave channel at the 0. 
H. Hindsdale Wave Research Facility (OHH-WRF) at Oregon State University. In 
order to determine the effects of wave groups on rubble mound structures, sea states 
with identical amplitude spectra but different wave group characteristics were simulated. 
The resolution of incident and reflected wave time series from the random wave 
simulation is required in order to evaluate the effects of wave groups. In order to 
quantify the two wave group characteristics proposed by Mase and Iwagaki (1986), the 
incident time series must be known. Medina et al. (1994) used a method developed by 
Kimura (1985) and Fassardi (1993) that extends an algorithm proposed by Goda and 
Suzuki (1976) to resolve the incident and reflected wave time series from closely spaced 
wave gauges. 



2 



CHAPTER I: Correlation of Wave Groups with Armor Damage 



Carstens, et al. (1967), first observed the effect of wave groups on breakwater 
damage when he analyzed two time series with significantly different groupiness, but 
similar spectra and wave statistics. The time series with the more grouped waves, taken 
from a field site in the Barents Sea, created more damage than the time series generated 
from the theoretical Neumann-spectrum. Johnson, et al. (1978), as a result of their 
laboratory studies on the effects of wave grouping on breakwater stability, concluded that 
wave grouping was an essential parameter in the model testing of rubble mound 
breakwater stability. Wave trains with similar spectra and wave statistics, but different 
groupiness, caused significantly different levels of damage, the greater damage being 
caused by the grouped wave trains. Burcharth (1979) also concluded that wave grouping, 
or the succession of the waves, is important in model design. 

Since the energy spectrum does not contain the phases of the individual wave 
components in the sea surface it is necessary to also resolve the phases of the individual 
wave components to reproduce an accurate succession of wave heights (Rye 1982). In 
his summary of breakwater design, Bruun (1989) recorded that the resolution of the 
phase spectrum identified apparent order where none is assumed in the Gaussian 
random wave model. He concluded that these ordered wave groups are more damaging 
than other waves. 

In contrast, van der Meer (1988) found that the spectral shape and grouping of 
waves had little or no influence on the stability of breakwater armor layers provided that 
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the average spectral period was used to simulate the sea state rather than the peak 
period. Hall (1994), in a study on bermed breakwaters, concluded that wave groups are 
only important during the stage of incipient motion, and that once motion is initiated 
wave groupiness is no longer important. 

Medina et al. (1990) found that much of the controversy could be resolved if the 
two independent parameters of wave groups proposed by Mase and Iwagaki (1986) were 
considered. One parameter measures the run length or the number of high waves in a 
time series, and a second parameter quantifies the magnitude of the variation of wave 
energy. In their laboratory tests, Medina et al. (1990) tested four different wave height 
time series proposed by Mase and Iwagaki (1986) shown in Figure 1.1, while keeping 
other secondary factors that affect breakwater stability constant. 




Figure 1.1 - Classification of wave height time series (Fassardi 1993) 



4 



Figure 1.1 shows that identical spectra may have different magnitudes of wave height 
variability, which may cause different levels of damage. Likewise, different spectra may 
have the same magnitude of wave height variability, and thus cause the same level of 
damage. Medina et al. (1990) simulated truncated Goda-JONSWAP spectra with 
different run lengths and different magnitudes of variation of wave energy according to 
Appendix B. 

Medina et al. (1990) concluded that the spectral shape and run length alone were 
not sufficient parameters to accurately describe wave groupiness, but that another 
parameter that measured the variation of wave energy more completely described the 
wave group. Medina et al. (1990) correlated an envelope exceedance coefficient a that 
measures the variation of wave energy with stability, and demonstrated that random 
wave trains with the same spectral shape may produce different levels of damage. 

Medina et al. (1994) found no significant correlation between the spectral peakedness y 
and armor damage, but concluded that the envelope exceedance coefficient a was 
necessary to evaluate the effects of random waves on structures in both physical and 
prototypical models. 
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CHAPTER II: Derivation of Envelope Exceedance Coefficient a 



From their analysis of storm waves, Mase and Iwagaki (1986) recommended using 
at least two parameters to quantify wave groups, one parameter that measures the 
magnitude of a sequence of high waves in a time series, and a second parameter that 
measures the magnitude of the variation of wave energy. Both parameters are shown 
schematically in Figure 2.1, where A(x,t) is the wave envelope function from the Hilbert 
transform defined by Medina and Hudspeth (1990). 




Figure 2.1 - Wave group parameters 
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The peak enhancement factor 7 of the Goda-JONSWAP spectrum (Goda 1985) 
controls the sharpness of the spectral peak and was selected to characterize the 
magnitude of the sequence of high waves in the time series. Rye (1980) related the 
magnitude of the sequence of high waves, run length, to the peak enhancement factor 7 
of the JONSWAP spectrum. The coefficient a was selected to characterize the 
magnitude of variation of wave energy. The envelope exceedance coefficient a is a 
dimensionless measure of the wave energy exceeding that associated with a prescribed 
wave height; e. g.,H mo . 

In a two dimensional wave channel, the wave height function H(x,t) is defined by 

H(x,t)= 2A(x,t) ( 2J ) 

If At = a constant sampling time interval, the discrete wave height function at a fixed 
location may be expressed as H(nAt). Defining a normalized measure of the variation of 
wave height as 

A H = H(nAt) - 1 (2.2) 

n tj ' 

n mo 

and the Heaviside step function (Farlow 1982) of AH n as 



U(AH n ) 



0 AH n <0 

1 ah" > 0 

n 



(2.3) 



then the magnitude of the variation of wave energy above H 1/10 of the time series is given 
as 
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(2.4) 



N 



«'4E (AH o ) 2 U(AH o ) 



N n=l 



where N is the number of discrete data points in the time series. Normalizing Eq. (2.4) 
by the expected value of a' yields the envelope exceedance coefficient a defined as 



a = 



a' 

EM 



(2.5) 



where E[cf] is defined as 



E[a']= | 




H 

TJ 

n i/io 




p(H) dH 



(2.6) 



where H is the wave height function at a fixed location H(t), and p(E[) is the Rayleigh 
probability density function of H(t) (Appendix A) defined by Chakrabarti (1987) as 



P(H) = 



H 

4m, 



■ exp 



- H 
8m A 



(2.7) 



and m 0 is the variance of the time series. If the following substitution is made (Sarpkaya 
and Isaacson 1981), 

H„ 10 = 5.091 < 2 ' 8 

Eq. (2.6) may be solved numerically to obtain Efc/] =0.001346. 
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CHAPTER ID: Resolution of Incident and Reflected Wave Time Series 



In order to determine experimental wave heights and to calculate envelope 
exceedance coefficients, the incident time series must be known. Because of reflection 
from models in a 2-dimensional wave channel, a method to resolve the incident and 
reflected random wave time series from wave gauges must be used. Bad data points in 
wave records caused by wave gauge malfunction are smoothed using the FORTRAN 
program CLEANREC that employs a smoothing spline package developed by H. J. 
Woltring (1986). A Fast Fourier Transform (FFT) analysis is used to compute the raw 
wave spectrum from which the incident and reflected random wave time series may be 
determined. Goda and Suzuki (1976) modified a method developed by Thornton and 
Calhoun (1972) to resolve the incident and reflected wave spectra from wave gauges 
spatially separated by a distance Af (see Fig. 3.1). 




O O 



y 






Figure 3.1 - Definition Sketch for two wave gauges 
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Kimura (1985) and Fassardi (1993) extended the Goda and Suzuki algorithm to 
include wave phases so that incident and reflected time series could be resolved. 

Goda and Suzuki (1976) analyzed two simultaneous wave records recorded by 
closely spaced wave gauges aligned in the direction of wave propagation by an FFT. If 
the incident t] and reflected f random wave time series at the i* wave gauge at frequency 
co are given by 





77 .= a cos(kx ; - cot+ £) 
f.= b cos(kXj+ ut+ j3) 


(3.1a) 

(3.1b) 


then the wave profile at the i* wave gauge may be expressed as 

(V r ; ) (x ^= AjCOso)t+ BjSinwt 


(3.2) 


where 


A,= a cos ^,+ b cos^ R 


(3.3a) 




B,= a sim/',- b sin^ R 


(3.3b) 




A 2 = a cos(kAf + i/',) + b cos(kAf + ^ R ) 


(3.3c) 




B 2 = a sin(kA£ + ^,)- b sin(kA£ + ^ R ) 


(3.3d) 


where 


the spatial phases of the incident and reflected waves, respectively, 

\p l = kXj+ £ 


are defined as 
(3.4a) 




^ R = kx,+ /3 


(3.4b) 




x 2 = x, + Af 


(3.4c) 


and X! 


and x 2 are the wave gauge positions as shown in Figure 3.1. The coefficients Aj 

10 



and Bj may be related to the complex-valued FFT coefficient X by (Dean and Dalrymple 
1984) 

X= Ai ~ jBi (3-5) 

2 

where j= \J - 1 . Goda and Suzuki (1976) solve Eqs. (3.3) simultaneously for estimates 

of the amplitudes a and b listed in Table 3.1. 

For irregular waves, the time series are the superposition of many wave 
components given by Eqs. (3.1). Because Eqs. (3.3) apply to each Fourier component in 
an irregular wave record, linear superposition may be used to resolve time series from a 
spectrum of Fourier components (Goda 1985). Goda and Suzuki (1976) observed that 
spectral estimates diverged (toward infinity) near frequencies satisfying the condition 
kAf =nir for n=0,l,2,..., because the term | sin kAf | in the denominator of the 
equations for the amplitudes a and b in Table 3.1 becomes small and errors from noise 
are amplified. Therefore, the wave gauge spacing Af determined the frequency limits of 
a band pass filter from which the waves could be separated into incident and reflected 
time series. The wave amplitudes a and b could be resolved effectively for the interval 
0. l7r<kAf <0.9x. Goda and Suzuki (1976) recommended the following effective band 
pass limits for experimental conditions: 

0.05L < Af < 0.45L . (3.6a) 

max min 

where L max and L min denote the wavelengths corresponding to the lower (f^) and upper 
(fm»x) frequency limits, respectively, of the band pass filter. Goda and Suzuki (1976) note 
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that the effective range may be taken slightly wider than given in Eq. (3.6a) such that 

0.03L < Af < 0.45L . (3.6b) 

max min 

Goda and Suzuki (1976) also recommended that the wave gauges be located at least one 
wavelength away from both the breakwater toe and the wave generator. 

Kimura (1985) extended the Goda and Suzuki algorithm to include situations 
where each incident wave component in the Fourier spectrum results in a different 
reflection coefficient from the structure. Kimura (1985) defined incident and reflected 
wave profiles by Eqs. (3.1), and the spatial phases by Eqs. (3.4a,b). Eqs. (3.1) are similar 
to Kimura’s Eqs. (9 &. 12) on pages 62-63 (1985) if the sign of the temporal term is 
changed from minus to plus in his Eq. (12). Kimura (1985) also solved for the spatial 
wave phases defined in Eqs. (3.4a,b). Values for a, b, \p u and i/'r may be calculated using 
the equations listed in Table 3.1 for Kimura. To improve the resolution of incident and 
reflected waves, Kimura recommended a slightly more conservative data window of 

0.15L < Af < 0.35L . (3-7) 

max min 

Fassardi (1993) followed the approach of Kimura (1985) to extend the Goda and 
Suzuki algorithm (1976). However, he defined the incident r? and reflected f time series 
at the i m wave gauge as 



77.= a cos(kx ; - ait — e) 


(3.8a) 


fj= b cos(kX;+ wt+ / 3 ) 


(3.8b) 



If the spatial phase of the incident and reflected waves, respectively, are defined as 

!/<,= kXj- £ (3.9a) 
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kx, + e 



(3.9b) 



rather than by Eqs. (3.4a,b), then the amplitudes and spatial phases are given by the 
equations listed in Table 3.1 for Fassardi (1993). The only difference between the 
Fassardi and Kimura algorithm is the sign of s in the definition of the incident time 
series phase. This sign change results in corresponding sign changes in the Fourier 
coefficients that are now given by 



A,= a cosi/q+ b cos\p K 


(3.10a) 


B, = a sini£j+ b sin^ R 


(3.10b) 


A,= a cos(-kAf+ ^,)+ b cos(kAf + \p K ) 


(3.10c) 


B, = a sin(-kAf + ^,)+ b sin(kAf + \^ R ) 


(3. lOd) 



These four equations are then solved simultaneously to obtain the equations for a, b, \p u 
and \p K listed in Table 3.1 for Fassardi (1993). 

In the experiments documented by Medina et al. (1994), the Fassardi algorithm 
was used to resolve the incident and reflected time series. Three sonic wave gauges 
aligned in the direction of wave propagation were centered 10m from the toe of the 
breakwater and spatially separated by A t = 1.22m according to Figure 3.2. 

Clean time series from each gauge were used in the FORTRAN program 
RECANSYG to resolve the incident and reflected time series. RECANSYG outputs the 
raw time series and spectra resolved from the wave gauge system in Figure 3.2. The 
program computes the time series mean and variance. As well, it computes the 
reflection coefficient, significant wave height, and mean period from the moments. The 
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Goda 






and 


a. = 


1 


2 I 


sinkAf | 


Suzuki 


b = 


1 


(1976) 


2 i 


| sinkAf | 


Kimura 


a = 


1 


(1985) 


2 1 


sinkAf | 



2 | sinkAf | 



\J(A 2 -AjCOskAf -B,sinkA f) 2 + (B 2 +A 1 sinkA£ -BjCOskA f ) 2 
\J( A 2 -A,coskAf +B 1 sinkAf) 2 + (B 2 -A,sinkAf -BjCOskAf) 2 



\J(A 2 ~A l coskAi -BjSinkA f) 2 + (B 2 +A,sinkAf -BjCoskA f) 2 
^(A^ -A,coskAf +B 1 sinkAf) 2 + (B 2 -A 1 sinkAf -BjCOskAf) 2 



\p= tan " 1 
\J/ K = tan' 1 



- A 2 + AjCoskAf + BjSinkAf 
B 2 + AjSinkAf - B,coskAf 

- A 2 + AjCOskAi - B r sirikA f 

- B 2 + A,sinkAf + B,coskAf 



Fassardi 

(1993) 



a= J( A,-A.coskAf +B.sinkA f) 2 + (B, -A. sinkAf -B.coskAf) 2 

2 | sinkAf | V 2 1 1 2 1 1 



b= i/(A- -A.coskAf -B.sinkA f) 2 + (B, +A. sinkAf -B.coskAf) 2 

2 | sinkAf | V 2 ' 1 V 2 1 1 



\p x - tan 
0 R = tan 



- A 2 + AjCoskAf - B,sinkA f 
B 2 - AjSinkAf- B,coskAf 

- A 2 + A,coskAf + B,sinkAf 
B 2 + A^inkAf- B,coskAf 



Table 3.1 - Summary of incident and reflected wave amplitudes and phases 



14 



incident and reflected time series were resolved using Gauges 1 and 2 and Gauges 2 and 
3 so that the algorithm could be verified by duplication. A sample run is analyzed. As 
explained in APPENDIX B - Simulation of Sea States, Run E1P1L7 represents a sea 
state simulated with the following parameters: 

El = Envelope 1 ; 7 = 10, a=2.52 
PI = Phase 1 ; 0 = 27 t/ 3 
L7 = Height 7 ; H s(7) =0.723 



AC = 1.22m 




Figure 3.2 - Sonic Wave Gauge Locations (Fassardi 1993) 



For a small data sample of ten seconds, Figure 3.3 shows how the clean time 
series are phase shifted at each wave gauge station. Figure 3.4 shows two minutes of the 
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clean time series at Gauges 1 and 2 and the resulting incident and reflected time series 
resolved by the Fassardi (1993) algorithm. The time series at Gauge 1 is identical to 
that at Gauge 2, except for the phase shift. The incident and reflected time series 
resolved from Gauges 1 and 2 are plotted on the same scale to show the relative 
magnitude of the incident and reflected wave amplitudes. The magnitude of the wave 
amplitudes of the reflected time series are about 25% of the incident amplitudes. Figure 
3.5 shows the composite spectrum at Gauge 2 and the incident and reflected spectra 
resolved from Gauges 1 and 2 and Gauges 2 and 3, respectively. The spectra in Figure 
3.5 have been smoothed using the non-statistical technique Box Car Averaging. The 
resolved spectra show that the energy in the reflected spectrum is about 5-10% of that in 
the incident spectrum. Breakwater reflection of about 25% observed in the time series 
results in the reflected spectral density being 5-10% of the incident spectral density, for 
each frequency component. The low frequency spike in the spectral density at 
approximately 0.04 Hz is attributed to a seiche in the wave channel caused by the start- 
up of the wave generator. The resolved spectra in Figure 3.5 are similar, indicating that 
the Fassardi (1993) algorithm produces the same spectra from different wave gauge 
pairs. At the truncation limits of 0.2140 Hz and 0.7644 Hz, the spectral density of the 
incident and reflected waves begins to diverge because the term | sin kAf | in the 
denominator of the equation for the amplitudes a and b in the Fassardi algorithm 
becomes small. Further examples of the resolution of incident and reflected time series 
and spectra from the output of RECANSYG for the four envelopes discussed in 
APPENDIX B are included in APPENDIX C - Time Series and Spectra. 



16 



Time Series - Run E1P1L7 




- - GAUGE 1 — GAUGE 2 - - - GAUGE 3 



Figure 3.3 - Phase lag between time series at three wave gauges 
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Figure 3.4 - Time series resolved from Gauge 1 and 2 
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Figure 3.5 - Composite and Resolved Spectra 
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SUMMARY AND CONCLUSIONS 



Multiple failures of existing rubble mound breakwaters designed with the Shore 
Protection Manual’s (1984) methodology have indicated a need to improve the design 
procedures for rubble mound structures. Many structural and environmental 
characteristics such as armor grading, structure permeability, wave period, storm 
duration, spectral shape, and wave groupiness have secondary effects on breakwater 
stability. Irregular waves generated solely from an energy spectrum with random phases 
do not accurately model real wave trains, and the succession of wave heights or 
groupiness is an important factor in the design of rubble mound breakwaters. The 
magnitude of the energy flux given by the envelope exceedance coefficient a is one way 
of quantifying the groupiness of a wave time series. Intermediate scale wave channel 
studies at Oregon State University documented by Medina et al. (1990), Fassardi (1993), 
and Medina et al. (1994) have demonstrated that the envelope exceedance coefficient a 
may be used to correlate wave groupiness with breakwater damage. No significant 
correlation was found between the peak enhancement factor 7 and rubble mound 
breakwater stability. Furthermore, the constant phase shift <£ applied to the time series 
for each wave envelope A(x,t) did not effect the armor damage. 

To calculate the envelope exceedance coefficient a from the incident wave time 
series, either the Kimura (1985) or the Fassardi (1993) algorithm may be used to resolve 
the incident and reflected waves from closely spaced wave gauges. The Fassardi (1993) 
algorithm, used by Medina et al. (1990) in their experiments at Oregon State University 
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is different from the Kimura (1985) algorithm by the sign of the random phase angle of 
the incident time series (Eq. 3.1a & 3.8a). This sign change does not have any effect on 
the value of a or the interpretation of the results. For the frequency range given by 
Table B.l.the Fassardi (1993) algorithm effectively resolves incident and reflected wave 
time series from closely spaced wave gauges aligned in the direction of wave 
propagation. Some divergence of the incident and reflected wave spectra is observed 
near the cut off frequencies of the truncated spectra. 

Since the wave resolution algorithms presented only resolve wave components in 
the frequency range determined by the wave gauge spacing, the spectra used to simulate 
the time series for the laboratory experiments were truncated accordingly. The variance 
of the full spectra was preserved for the truncated spectra by considering the squared 
significant wave height ratio given by Eq. (B.6). For the physical simulations performed 
by Medina et al. (1990), phase spectra resulting in the highest and lowest values for the 
envelope exceedance coefficient a w r ere chosen for each amplitude spectra. Significant 
wave heights of consecutive runs were increased in discrete increments such that the 
stability numbers for the two armor rock sizes were equal in consecutive runs. 
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APPENDIX A: Rayleigh Distribution 



The Rayleigh probability density function (pdf) for the wave amplitudes a ( = 
H/2) is given by 



P (a) = —exp 
“r 



2 



a 

CC n 



U (a) 



a R >0 



(A. la) 



where or R = Rayleigh parameter (Hoffman and Karst 1975) and U(a) = Heaviside step 
function. The Rayleigh cumulative distribution function (cdf) for wave amplitudes is 



P(a) = 1 - exp 



1 

2 



a 



U (a) 



The following 4 forms of (A. la) may be found: 



P(a) = -^exp 



/ux 2H 
P(H) = -7-exp 
H 



P( a ) = — exp 
m 



, YJ V H 

P (H ) = - — exp 
4 m 

0 



- 




- 






2 


t r 


a 




— 







4 


p. 









2 




H 






H 

rms 





2 



1 


H ] 


2 


2 







U (a) 



U (H ) 



a R >0 



U(a) ; ^>0 



U(H) ; H >0 



; m o >0 



m >0 

o 



(A. lb) 



(A. 2a) 



(A. 2b) 



(A. 2c) 



(A. 2d) 



where = average of the amplitude a; H ras = root-mean-square wave height; and m 0 = 
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variance of the time series for the water surface elevation 7j(t). When the random 



variable a in (A. 2) is scaled by the standard deviation of the time series 




, then the 



following change of variables gives: 




p«) - P(a)4r 
d£ 



Ptt) = 




P(0 = £exp 



e 

2 



u«) 



(A.2e) 



Eqs. (A. 2) may be derived from (A. la) by solving for the generic Rayleigh parameter a R . 

The percent or fraction of wave amplitudes greater than wave amplitude 
a,, (= H 0 /2), may be computed from (A. la) by 



oo 

n = | p(a)da 



-> 2 




(A. 3) 



and the natural logarithm of (A. 3) gives 
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il = i/-21n (n) 



(A.4) 



o?„ 



where 0<n<l. The average wave amplitude greater than a n may be determined from 

co co 

a n | p (a)da = j ap(a)da 

a t 

that may be integrated by parts to obtain 



a „(n) = na RV /-21n(n) + | exp 
by (A. 3) and (A.4). Finally, integration gives 



2 

2 






da 



a , i 

-2. = v/-21n(n) + - 






— l-erf( 
2 L v 


v/-ln(n)|J 


(A. 5a) 




\ 


(A. 5b) 


1 erfc(v/- 


In (n) ) 



= i/-21n(n) + I 
n 

where erf( • ) = the error function and erfc( • ) = complementary error function 

(Kreyszig 1983). In order to relate the Rayleigh parameter a R to the average pi,, set n = 
1 in Eq. (A. 5) and 



a „ a . M 

n _ 1 * a 



a R “r a R 



X 

2 



a R = M. 

and substitution into (A. la) gives (A. 2a). 



(A. 6) 



2 

x 



24 



Alternatively, (A. 5) may be derived directly from special functions according to 



a n | p(a)da = j ap(a)da 



= 



- - v/2 r 

2 



a 

-, 0 ,— 

2 2<* 2 



where the generalized incomplete Gamma function r(a,z 0 ,z 1 ) (Wolfram 1991) is defined 



as 



r(a,z 0 ,z ] ) = r(a,z 0 ) - r(a,z,) 

and the incomplete Gamma function T(a,z) (Wolfram 1991) is defined as 

CD 

f(a,z) = f(a) - G(a,z) = j t ,_, exp( -t) dt 



G (a,z) = 



t *''exp(-t) dt 



Substituting for r(a,z 0 ,z,) 






2 ’ 2 al 



_£ 

2 



-rfi 






i/T exp( -t) dt 



= f v/t~exp( -t) dt 
^ 1 






= n/-ln(n) + 1 -erf(v/-ln(n) ) 



= n i/ — In ( n ) + -^-erfc(^-ln(n) ) 
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giving (A. 5). 

To relate the Rayleigh parameter a R to H^, make the following change of 
variables in (A. la): 

y = a 2 ; dy = 2 a da = 2\/y"da 



and 

P (a) da = p (y)dy 



giving 



p(y) = p(a)~ 

dy 



exp 



2«l 



_ y 



2 a; 



U(y) 



(A.7) 



The average of a 2 may now be computed from 



’2 _ 



2 = 
'mss 



a_i_ = y 



yp(y)dy 






exp 



_ >' 



2« 2 r 



dy 






(A. 8) 



and the Rayleigh parameter a R may be defined by 
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H 



a R 



\Jl 2j2 




Substituting (A. 6) and (A. 9) into (A. 5b) gives 




— erfc(v/-ln(n) ) 
n 



H 

C 

W 



= \/-ln(n) + -^erfc(i/ -In (n) ) 
2n 



-?i- = v/-81n(n) + ^Lefrc( v Mn(n) ) 

/ n 

V m o 



The significant wave height H, or average of the highest 1/3 wave heights I 
computed from (A. 5b) for n = 1/3 (dropping the overbar average notation) 

= JL = 2.00215 

2or R 2a R 

2a - H >( H '«) 

R 2.00215 

Substituting (A. lib) into (A. 5b) gives 



(A. 9) 



(A. 10a) 

(A. 10b) 

(A. 10c) 

1/3 may be 
according to 

(A. 11a) 
(A. lib) 
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(A. 12) 



H 






= (2.00215) 



-ii 



i/-21n(n) + — 



j erfc(v/-ln(n) ) 



To find the mode (or most probable) value of a^,^ (= H^^/2), the maximum value of 
(A. la) occurs when 



dp(a) = n . 
- U > 
da 



a = a 



mode 



(A. 13) 



giving 



a mode ^ 



(A. 14a, b) 



and 



mode 



2 

T 



H 



mode 



H 



_ 1 _ 

& 



(A.15a,b) 



H 



mode 



= 2 ; 



H 



mode 






■= (2.00215) 



-1 



(A.15c,d) 



The mean or median of a occurs when P(a) = 0.5 in (A. lb). The natural logarithm of 
(A. lb) for P(a mc<1 ) =0.5 is 

— = i/21n (2) (A. 16) 

a R 
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giving 



a m«i 



41n (2) 

7T 



— = \/ In ( 2 ) 



(A.17a,b) 



— — = J 81n(2) 






'/2ln(2) 

2.0015 



(A. 17c, d) 



These results are summarized in Table A.l 



height 




H 


H 


H 


H 




H 

rms 






H . 

mode 


> 


=t 1 oo 


i 

If 


2.0 


(2.00215) 


H msd (P = l/2) 












4 

> 


ln(2) 

7T 


i/ln(2) 


v/81n(2) 


v/21n(2) 

2.00215 


H t ( = mean ) 


2.0 


& 

2 


/ 2tt 


\J v/ 2 

2.00215 




3.19497 


1.41573 


4.0043 


1.0 


TT 

^1/10 


4.06198 


1.79992 


5.09094 


1.27137 


i_r 

XX 1/100 


5.32423 


2.35924 


6.67293 


1.66644 



Table A.l -Rayleigh distributed wave height relationships 
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APPENDIX B: Simulation of Sea States 



In order to correctly describe wave groups in laboratory studies, the incident wave 
train must be known. An algorithm developed at Oregon State University by Fassardi 
(1993) was used to resolve the incident and reflected wave trains in a finite bandwidth 
determined by the wave gauge spacing Af (Table 3.1). One broad (7=1) and one 
narrow (7=10) truncated Goda-JONSWAP spectrum (Goda 1985) were simulated 
according to 






C.KLH*, 

1 R s.lrnnc 



f 4 f cxp[-d-l)*/2*‘] 

JLexp[-1.25(I)- 4 b '' 



(B.l) 



'f 



0.0624 

0.230+ 0.03367- 0.185(1.9+ 7)- 1 



0= 



0.07 : f<f 
0.09 : f>f 

P 



(B-2) 

(B.3) 



for f min < f < f max , where f p is the peak frequency, H R is the squared significant wave height 
ratio defined in Eq. (B.6), and H s lrunc is the significant wave height of the truncated 
spectrum. Eq. (B.l) is different than Eq. (30) in Fassardi (1993) in that Fassardi (1993) 
does not document the use of the constant c, defined by Eq. (B.2). The spectra were 
truncated according to 



f . 

min 



0.7f 

p 



f = 



max 



2 - 5f P 



(B.4) 
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To define the truncated spectra and preserve the variance of the full spectra, Eq. 
(B.l) includes the dimensionless coefficient H R . If the significant wave height is given in 



terms of the variance (Table A.l) by 



H =4.0043 

5 




(B.5) 



then the squared significant wave height ratio H R may be defined as 




(B.6) 



where H s>fu „ is the significant wave height of the full spectra. For peak frequencies 
between 0.1 and 0.5 Hz, the average dimensionless variance preserving coefficients for 
7=1 and 7 = 10 are H R =1. 0375 and H R =1.0133 respectively. 

The mean frequency of the truncated spectra defined by the first and zeroeth 
spectral moments 

? = — (B.7 

m 0 

was held constant at ?= 0.3333 Hz such that the peak frequency of the truncated spectra 
could be computed from the frequency ratio 



f 

For peak frequencies between 0.1 and 0.5 Hz, the average frequency ratios were 




(B.8) 
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f R =0.81 15 and f R =0. 9173 for 7 = 1 and 7= 10, respectively. Accordingly, the peak 
frequencies used to simulate the truncated spectra from Eq. (B.l) were f p =0.2705Hz 
and f p = 0. 3058Hz for 7 = 1 and 7= 10, respectively. 

The wave gauge spacing defined by Figure 3.1 was 1.22m, and the water depth 
was 3.05m (Fassardi 1993). If the band pass filter limits proposed by Goda and Suzuki 
(1976) given by Eq. (3.6a) are used, and the linear wave theory dispersion relationship is 
given by 



CJ = 



2tt , 27rh 
g — tanh 

L L 



(B.9) 



the minimum and maximum resolvable frequencies for A? = 1.22m are f mIn =0. 1299 Hz 
and f m „=0.7589 Hz. For cases where 7 = 1 , the minimum and maximum frequencies of 
the truncated spectra are f min =0. 1894 Hz and f m „ = 0.6762 Hz. For 7 = 10, the minimum 
and maximum frequencies of the truncated spectra are f min =0.2140Hz and f m „=0.7644 
Hz. These frequency limits for the truncated spectra are reasonably close to those 
resolvable by wave gauges spatially separated by A£ = 1.22m. 

Different phase spectra produce time series with different wave grouping 
characteristics. Accordingly, time series were synthesized using two different phase 
spectra for each of the truncated Goda-JONSWAP spectra given by Eq. (B.l) with 

7 =l,l 0 and ?=0.3333 Hz. The phase spectra were chosen from 100 random 

Deterministic Spectral Amplitude (DSA) simulations for each amplitude spectra. Values 
of a were calculated from each DSA simulated time series according to Eq. (2.5). For 
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each 7 , the phase spectra that produced the wave envelope A(x,t) with the highest and 
lowest wave height variability characterized by a were chosen for the physical 
simulations. The parameters used for wave simulation are summarized in Table B.l. 



7 


Cl 


H r 


f R 


fp 


resolvable for : 
A? = 1 . 22 m 
from Eq. (3.6a) 


f . 

rrun 


f 

A max 


Low a 
from 
DSA 


High a: 
from 
DSA 






































1 


0.3123 


1.0375 


0.8115 


0.2705 


0.1299 


0.7589 


0.1894 


0.6762 


0.51 


2.18 


10 


0.1134 


1.0133 


0.9173 


0.3058 


0.1299 


0.7589 


0.2140 


0.7644 


0.23 


2.52 



Table B.l -Summary of simulation parameters for At = 1.22m and h=3.05m 

(frequencies in Hz) 



Realizations with the same envelope, but different phase-shifted wave profiles are 
given by (Medina et al. 1994) 



M 

^(x,t) = £a m cos[kx-«t-(e+$)] 

m = 1 



(B- 10 ) 



where <f> is a constant phase shift applied to all M wave components. Three different 
phase shifted realizations (<£=0,27 t/ 3, 4rr/3) and the replicate (<fr=2ir) for each of the 
four wave envelopes were tested to determine if realizations shifted by a constant phase 
produce different levels of damage. For each y, A(x,t), and constant phase shift 4>, a 
sequence of time series were tested with increasing significant wave heights. 
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The breakwater model tested by Medina et al. (1990) was divided in the direction 
of wave propagation into two sides with a different armor layer rock size (viz., W L = 128.5 
N and W s =99.1 N) on each side. The monochromatic wave height corresponding to the 
zero damage condition (SPM 1984) can be calculated from 

i 

3 

(B. 11) 





r 


3 


wk d 


^1-1 

Pw 


cot9 


p r 



where W is the median value of the mass distribution of rocks in the armor; K D is the 
stability coefficient; 0 is the angle of the breakwater slope measured from horizontal, 
and p r and p w are the weight densities of the rocks and water, respectively. Solving Eq. 
(B.ll) numerically for the small rock size (W s =99.1 N) yields a monochromatic wave 
height of H = 0.55m for the zero damage condition. Breakwater model and armor rock 
characteristics are given in Table B.2. 



W L (N) 


W s (N) 


Kd 


p f (kN/m 3 ) 


Pw (kN/m 3 ) 


cot9 


128.5 


99.1 


4 


27.4 


9.8 


2 



Table B.2 - Summary of breakwater model characteristics (Fassardi 1993) 



For the zero damage condition, if the representative design wave height for a 
random sea is H 1/10 , then the design significant wave height for the small rock size is 
determined (Table A.l) by 
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(B- 12) 



H s =- Q,55 ^-= 0.43m for W s = 99.1 N 



5 1.27137 



An empirically derived number which controls the stability of the armor layer (Medina et 



al. 1990) is given by 



N = 

s 





3 



(B. 13) 



where N s is the stability number, and H s is the design significant wave height. The 
significant wave heights in the physical simulations were increased in discrete increments 
such that the stability numbers given by Eq. (B.13) for the large (WJ and small (W s ) 
rocks were equal in consecutive runs. Seven significant wave heights were tested for 
each of the sixteen realizations. The significant wave heights varied from the design 
significant wave height of the small rock size to the maximum wave height, avoiding 
breaking. From Eq. (B.13), the significant wave height of the k* run is given by 



Thus, twenty-eight realizations of four phases and seven wave heights were 
simulated for each of four target envelopes shown in Figure B.l. For each realization, 
runs of approximately 30 minutes long containing N=2 15 =32768 points sampled at 
At=0. 06 seconds were simulated. 



r 



W, -r 



(B. 14) 



H . =0.43 

s.(V) 
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Figure B.l - Simulation of Realizations with Wave Height H. 



APPENDIX C: Wave Time Series and Spectra 



Section Page 

1. Run E1P1L7 38 

2. Run E2P1L7 40 

3. Run E3P1L7 45 

4. Run E4P1L7 50 
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RESOLVING INCIDENT AND REFLECTED TRAINS FOR RECORD: elp!17 



Gauge No. : 1 

Raw Time Series Mean Value : 
Gauge No. : 2 

Raw Time Series Mean Value : 
Gauge No. : 3 

Raw Time Series Mean Value : 



0.009327 


Variance : 


0.040746 


0.010709 


Variance : 


0.038375 


0.006861 


Variance : 


0.036691 



RESOLVED INCIDENT AND REFLECTED WAVES 



Analysis from Gauges 2 and 3 

Incident wave Spectrum’s Variance: 0.032890 

Reflected wave Spectrum’s Variance: 0.001678 

Breakwater’s reflection coefficient : 0.2259 
Incident significant wave height : 0.7261 
Reflected significant wave height : 0.1640 
Mean Period (T01) from incident spectrum : 3.065 

Incident wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.032890 

Reflected wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.001678 

Analysis from Gauges 1 and 2 

Incident wave Spectrum’s Variance: 0.033244 

Reflected wave Spectrum’s Variance: 0.002225 

Breakwater’s reflection coefficient : 0.2587 
Incident significant wave height : 0.7300 
Reflected significant wave height : 0.1889 
Mean Period (T01) from incident spectrum : 3.074 

Incident wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.033244 

Reflected wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.002225 
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Resolved from Gauges 2,3 (m) Resolved from Gauges 2,3 (m) Gaige 3 (m) Gauge 2 (m) 



Timo Sei-icss - 3R-vi n E1P1L7 




Time: C s ) 



Time Series — R.un E1P1L7 




O 30 <50 S>0 120 

Time C s ) 



Incident T ime Series - Run E1P1L7 




O 30 <50 S>0 120 

Timer C s ) 



ReHectecI Time Series - R_ti Ell D 1L7 




Time C s ^> 
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RESOLVING INCIDENT AND REFLECTED TRAINS FOR RECORD: e2p!17 



Gauge No. : 1 

Raw Time Series Mean Value : 0.009919 Variance : 

Gauge No. : 2 

Raw Time Series Mean Value : 0.003559 Variance : 

Gauge No. : 3 

Raw Time Series Mean Value : -0.010094 Variance : 

RESOLVED INCIDENT AND REFLECTED WAVES 

Analysis from Gauges 2 and 3 

Incident wave Spectrum’s Variance: 0.031342 

Reflected wave Spectrum’s Variance: 0.002484 

Breakwater’s reflection coefficient : 0.2815 
Incident significant wave height : 0.7089 
Reflected significant wave height : 0.1996 
Mean Period (T01) from incident spectrum : 3.245 

Incident wave train, Mean and Variance from TS(2,3) 

TS Mean : 0.000000 TS Variance : 0.031342 

Reflected wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.002484 

Analysis from Gauges 1 and 2 

Incident wave Spectrum’s Variance: 0.033780 

Reflected wave Spectrum’s Variance: 0.002935 

Breakwater’s reflection coefficient : 0.2948 
Incident significant wave height : 0.7359 
Reflected significant wave height : 0.2169 
Mean Period (T01) from incident spectrum : 3.244 

Incident wave train, Mean and Variance from TS(1,2) 

TS Mean : 0.000000 TS Variance : 0.033780 

Reflected wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.002935 



0.039530 

0.041591 

0.041708 
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Time Series - TS-vm E2P1L7 




T i m e C s ) 



Time Series — ITun E2P1L7 






13 

"o 

ig 




Time: <» 



Incident Time Series — ITun E2P1 ^ 




Time: 



ITe llected Time Series — ITu n E2P1 LV 
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oo 



1 § 

“o 

5£3 

Dei 




Time C s ) 
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Time Series — K.un E2P1 LT 



Op 




Time: ( s ) 

Time Series — !Rnn E2P1L7 




Time C s ^> 

Incident Time Series — Run E2P1 LT 







Time 

Rellectecl Time Series — Run E2P1L7 
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Time C s ) 
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Resolved from Gauges 2,3 (nf 2/Hz) Resolved from Gauges 1 ,2 (nft/Hz) 



Composite Spectrum - Run E2P1L7 




Composite Spectrum — Run E2P1 L7 



go 

• 7 ^ 







Inc. nncl Ref. Spectra — Run E2P1 E7 




rN T CEDHN*T 



Inc. mul Ref. Spectra — Ru n E2P1 L7 




EN'CEDENT 



IOmLECTEX> 
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Time Series - Run E2P1L7 




- - GAUGE 1 — GAUGE 2 - - - GAUGE 3 
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RESOLVING INCIDENT AND REFLECTED TRAINS FOR RECORD: e3p!17 



Gauge No. : 1 

Raw Time Series Mean Value : 0.011743 Variance : 

Gauge No. : 2 

Raw Time Series Mean Value : 0.002714 Variance : 

Gauge No. : 3 

Raw Time Series Mean Value : -0.002450 Variance : 

RESOLVED INCIDENT AND REFLECTED WAVES 

Analysis from Gauges 2 and 3 

Incident wave Spectrum’s Variance: 0.035032 

Reflected wave Spectrum’s Variance: 0.001793 

Breakwater’s reflection coefficient : 0.2262 
Incident significant wave height : 0.7494 
Reflected significant wave height : 0.1695 
Mean Period (T01) from incident spectrum : 3.056 

Incident wave train, Mean and Variance from TS(2,3) 

TS Mean : 0.000000 TS Variance : 0.035032 

Reflected wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.001793 

Analysis from Gauges 1 and 2 

Incident wave Spectrum’s Variance: 0.034591 

Reflected wave Spectrum’s Variance: 0.002552 

Breakwater’s reflection coefficient : 0.2716 
Incident significant wave height : 0.7447 
Reflected significant wave height : 0.2023 
Mean Period (T01) from incident spectrum : 3.078 

Incident wave train, Mean and Variance from TS(1,2) 

TS Mean : 0.000000 TS Variance : 0.034591 

Reflected wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.002552 



0.041173 

0.038172 

0.038333 
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Resolved from Gauges 1 ,2 (m) Resolved from Gauges 1,2 (m) Gauge 2 (m) 



Time Series - K.un E3P1 L7 




Time ( s ) 



Time Series - Run E3P1 L7 




Time ( s ) 



Incident Time Series — Run E3P1L7 
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Resolved from Gauges 2,3 (m) Resolved from Gauges 2,3 (m) Gouge 3 (m) Gauge 2 (m) 



Time Series — K-U n E3P1L7 
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Time Series — Run E3P1L7 




O 30 60 PO X 2.0 

Time C s ) 



Incident Time Series — Run E3P1R7 




O 30 <50 90 120 

Time C s ) 



Re llected Time Series — Run E3 PI L,7 



T 








* 








- 








A A JTi. . f\ A A A 




A A. -v — /I A AA 


r A f\ A « A A A. A A 








‘Vvv^Wvv \ 


- 
















- 






> < 



O 3 0 <5 0 90 120 

Time C s ) 



47 



Composite Spectrum - Run E3P1L7 




Composite Spectrum — Rxm E2P1L7 
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Time Series - Run E3P1L7 




Time (s) 



- - GAUGE 1 — GAUGE 2 - - - GAUGE 3 



49 



RESOLVING INCIDENT AND REFLECTED TRAINS FOR RECORD: e4p!17 



Gauge No. : 1 



Gauge No. : 2 



Gauge No. : 3 



0.015713 


Variance : 


0.043060 


0.025821 


Variance : 


0.041359 


0.029771 


Variance : 


0.044257 



RESOLVED INCIDENT AND REFLECTED WAVES 



Analysis from Gauges 2 and 3 

Incident wave Spectrum’s Variance: 0.036691 

Reflected wave Spectrum’s Variance: 0.002226 

Breakwater’s reflection coefficient : 0.2463 
Incident significant wave height : 0.7670 
Reflected significant wave height : 0.1889 
Mean Period (T01) from incident spectrum : 3.206 

Incident wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.036691 



Reflected wave train, Mean and Variance from TS(2,3) 
TS Mean : 0.000000 TS Variance : 0.002226 



Analysis from Gauges 1 and 2 

Incident wave Spectrum’s Variance: 0.036596 

Reflected wave Spectrum’s Variance: 0.002472 

Breakwater’s reflection coefficient : 0.2599 
Incident significant wave height : 0.7660 
Reflected significant wave height : 0.1991 
Mean Period (T01) from incident spectrum : 3.212 

Incident wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.036595 

Reflected wave train, Mean and Variance from TS(1,2) 
TS Mean : 0.000000 TS Variance : 0.002472 
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Resolved from Gauges 1,2 (m) Resolved from Gauges 1,2 (m) Gauge2(m) Guugc 1 (m) 



Tixiies Series — Run E4P1 L7 




Time C s ) 
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Resolved from Gauges 2,3 (m) Resolved from Gauges 2,3 (m) Gauge 3 (m) Gauge 2 (m) 



Time Series — Run E4P1 L7 
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Resolved Rom Gauges 2,3 (nf2/Hz) Resolved from Gauges 1,2 (nfl/Hz) 



Composite Spectrum — Run E*4P1 JL7 




Composite Spectrum — Run E4P1L7 
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Time Series - Run E4P1L7 




- - GAUGE 1 — GAUGE 2 - - - GAUGE 3 
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APPENDIX E: Notation 



A(x,t) 

Aj 

^TTl&d 

^modc 



‘h-ms 

a 

B, 

b 

Cl 

da 

dH 

d£ 

E(a') 
erf( • ) 
erfc( • ) 
f 
? 



= wave envelope function 
= FFT coefficient at i* wave gauge 
= the median of a 
= most probable value of a 

= the average wave amplitude greater than a„ 

= root mean square wave amplitude 
= incident wave amplitude 
= FFT coefficient at i Ih wave gauge 
= reflected wave amplitude 

= Goda-JONSWAP spectrum parameter defined in Eq. (B.2) 
= wave amplitude interval 
= wave height interval 
= interval of £ defined in Eq. (A.2e) 

= expected value of a' 

= error function of ( • ) 

= complementary error function of ( • ) 

= wave frequency 

= mean frequency of spectrum 



f, 



nux(min) 



= maximum (minimum) cut-off frequency 
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f = peak frequency of spectrum 

f R = frequency ratio given by Eq. (B.8) 

G(a,z) = gamma function from APP. A 

o = gravitational constant 

H = wave height 

H(t) = wave height function at a fixed location 

H(x,t) = wave height function 

H(nAt) = discrete wave height function at x=0 

H mcd = the median of H 

H modc = the most probable value of H 

H r = squared significant wave height ratio given by Eq. (B.5) 

H,,,,, = root mean square wave height 

H s = significant wave height = H 1/3 

H sful | = significant wave height of full spectrum 

H s>(k) = significant wave height of the k* run 

H s>tninc = significant wave height of truncated spectrum 

H! = average of the wave heights 

H 1/3 = average of the highest one-third of the wave heights 

H 1/10 = average of the highest one-tenth of the wave heights 

H]/ioo = average of the highest one-hundredth of the wave heights 

h = water depth 

j = f- 1 
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K D 


= armor layer stability coefficient 


k (=2*/L) 


— wave number 


L 


= wave length 


L m ix(min) 


= maximum (minimum) wave length 


m 0 


= variance of the time series, zero-th spectral moment 


mi 


= first spectral moment 


n 


= percent or fraction of wave amplitudes greater than a„ 


N 


= number of data points 


N $ 


= stability number defined in Eq. (B. 12) 


P(-) 


= Rayleigh cumulative distribution function of ( • ) 


P(0 


= Rayleigh probability density function of ( • ) 


s,(f) 


= Goda-JONSWAP spectrum defined by Eq. (B.l) 


T 


= wave period 


t 


= time 


U( •) 


= Heaviside step function of ( • ) 


w 


= the median value of the mass distribution of rocks in the armor 


W L 


= large armor layer rock size 


W s 


= small armor layer rock size 


X 


= complex valued FFT coefficient 


Xj 


= position of i 0 * wave gauge from wave maker 


y 


= change of variable from Eq. (A. 7) 


a 


= envelope exceedance coefficient 
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a' 


= magnitude of the variation of wave energy above Hi /J0 


a K 


= Rayleigh parameter 


0 


= reflected wave phase 


r(a,z) 


= incomplete gamma function from APP. A 


r(a,Zo,z,) 


= generalized incomplete gamma function from APP. A 


7 


= peak enhancement factor 


AH n 


= normalized measure of the variation of wave height 


Af 


= wave gauge spacing 


At 


= sampling time interval 


£ 

f« 


= incident wave phase 
= reflected time series at i* wave gauge 


l?i 


= sea surface elevation 
= incident time series at i* wave gauge 


e 


= angle of the breakwater slope measured from horizontal 


Mi 


= average of the amplitude a 


$ 


= change of variable from Eq. (A.2e) 


Pr 


= weight density of rock 


Pw 


= weight density of water 


a 


= Goda-JONSWAP spectrum parameter defined in Eq. (B.3) 


<3? 


= constant phase shift applied to each realization 


'/'ICR) 


= incident (reflected) spatial wave phase 


u (=2 tt/T) 


= circular wave frequency 
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